Phase Diagram of Magnetization Reversal Processes in Nanorings 
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We investigate magnetization reversal processes of magnetic nanorings. Using a recently developed 
efficient cartesian coordinates fast multipole method, magnetization switching phase diagrams are 
constructed for such structures. We find that the dominant switching mechanisms strongly depend 
on the relative spatial dimensions of the rings, i.e. height, width and radius. Furthermore, they are 
found to depend on a characteristic exchange length scale, which is determined by the competing 
short-range ferromagnetic and long-range dipolar interactions in these systems. Detailed simulations 
allow us to identify three novel switching mechanisms, which should be useful for building high 
density storage systems. 



I. INTRODUCTION 

Mesoscopic magnetic elements have attracted tremen- 
dous attention during the last two decades .— ^ In partic- 
ular, there has been an upsurge of interest in small mag- 
netic rings during the last few years.— & These systems 
exhibit many novel physical phenomena with promise 
for applications, such as high density data storage^^, 
magnetic logics^ and biomedical sensorin g 10 ' 11 . Further- 
more, nanoring geometries have proven to be a fertile 
platform for the investigation of fundamental physical 
questions concerning domain walls, — One of the most 
useful properties of magnetic rings is that they can sta- 
bilize flux-closed vortex states. These vortex states are 
free of stray fields and largely insensitive to edge imper- 
fections, thus resulting in highly controllable magnetiza- 
tion reversal processes. Vortex states have been found to 
be stable in rings with diameters as small as 10 nm. In 
contrast, vortex states in magnetic discs exist only when 
their size is fairly large, i.e. with diameters larger than 
100 nm, due to the existence of vortex core regions, which 
lead to high energy penalties^ 

Experimentally, nanomagnctic rings have been pro- 
duced with a wide range of spatial dimensions, extend- 
ing from dozens of micrometers down to nanometer 
scales^. Various states have been observed to be sta- 
ble or metastable at remanence, i.e. in the absence of 
an applied magnetic field, including the vortex state (V), 
the onion state (O), i.e. two ferromagnetic domains sep- 
arated by domain walls at opposing ends of the ring 
structure, and the twisted (or saddle) state (T)i^), i.e. 
a vortex state interrupted by a 360° domain wall. Sev- 
eral works have been devoted to determining the phase 
diagram for nanoring structure o 13 ' 16 i 17 ' 18 . Apart from 
the phase diagram at remanance, magnetization rever- 
sal processes driven by external magnetic fields are an- 
other important characteristic property of magnetic el- 
ements. Several switching processes have already been 
discovered^. Experimentally, switching phase diagrams 
have been constructed for micron-scale rings^ However, 
the properties of nano-scale rings are still under investi- 
gation, and a number of new features have recently been 
discoveredi2ii22i22 

In this paper, we identify several switching processes, 



which according to our knowledge have not yet been re- 
ported in the literature, and construct switching phase 
diagrams for magnetic nanorings. We use Monte Carlo 
simulations combined with a recently developed efficient 
cartesian coordinates fast multipole tcchniqu o 24 i 25 i 26 i 27 
to analyze the magnetization reversal processes for vari- 
ous nanoring dimensions. Neglecting the crystalline en- 
ergy, the total energy (E) of a magnetic nanoparticle in 
a magnetic field consists of three terms: the exchange 
interaction, the dipolar interaction, and the Zeeman en- 
ergy, 



E = J £ &-£-+^£ 

<i,j> 
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Here J > is the ferromagnetic exchange constant, which 
is assumed to be non-zero only for nearest neighbors, D 
is the dipolar coupling parameter, and is the displace- 
ment vector between lattice sites i and j. Note that S is 
a dimensionless unit spin vector with magnetic moments 

Experimentally, the most studied magnetic nanoma- 
tcrials are made of Co, Fc, permalloy and supcrmalloy. 
For these materials, the ratio D/Ja 3 falls into the range 
between 10~ 3 and 10~ 5 , where a ~ 0.3 nm is the lat- 
tice constant. In the Monte Carlo simulations, single 
spin updates are adopted to yield quasi time dependent 
behavior^ We choose Ja 3 / D = 5000 which is close to 
the parameters of Co. Since we are trying to discover the 
generic property of nanomagnets, the absolute values of 
parameters used in this paper are not as important as 
their relative magnitudes. To state our conclusions in 
the most general terms, we define an exchange length, 
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whose overall prefactor a may vary between 0.5 and 2, 
depending on the material. Here we choose a = 1, and 
assume a typical value a = 0.3 nm, such that L ex ~ 20 
nm. This exchange length is very important in deter- 
mining the properties of nanomagnetic systems. It is a 
measure of the spatial extent of magnetic domains. For 
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instance, the typical diameter of vortex cores has been 
shown to be 1.3L ex r^- The temperature for the Monte 
Carlo simulations is set to 50K. and at least 10 4 Monte 
Carlo steps are performed for each field point. To facil- 
itate the following discussion, we define R as the outer 
ring radius, Ri as the inner ring radius, w = R — R t 
as the width and h as the ring height (see the inset of 
FigfJta)). 
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FIG. 1: (color online) Switching phase diagrams of magnetic 
nanorings. As indicated in the inset, h denotes the ring height, 
Ri and R are the inner and outer ring radii, and w = R — Ri 
is the width. A magnetic field H is applied in plane horizon- 
tally. The magnetic exchange length L ex is defined in the text. 
Phase diagrams are shown for the cases (a) 0.3L ex < h < L ex , 
(b) h > Lex, (c) h < 0.3L ex . The shaded areas in the regime 
of small w's represent rings whose remanence states are out- 
of-plane ferromagnets, and hence not relevant for switching. 
The color filled areas represent the different switching pro- 
cesses discussed in the text, green: onion - onion (O-O); red: 
onion - twisted - onion (O-T-O); yellow: vertical onion - vor- 
tex - onion (O v -V-0); gray: onion - twisted - vortex - onion 
(O-T-V-O); cyan: onion - vortex - vortex core - onion (O-V- 
VC-O); white: onion - vortex - onion (O-V-O). The point 'A' 
indicates the onion - vortex transition for discs. This critical 
point moves towards larger (w, R ) values when h decreases. 



netization reversal processes which have been intensively 
studied and discussed in the literature. 



A. Previously Reported Switching Processes 

1. The first process is one-step (or single) switching 
O-O, whose magnetization curve M(H) is illus- 
trated in FigfJ^b)). This is a direct onion-to-onion 
switching process, whereby an onion state is re- 
versed to the opposite onion state directly via a 
coherent rotation, as two magnetic domain walls 
move in the same rotational direction. It is worth 
noting that when the ring radius is sufficiently small 
(R < L ex ), the onion state is replaced by a single 
domain ferromagnetic state. 

2. The second process is two-step or double switch- 
ing O-V-O, shown in Fig. ^c). Here, two oppo- 
site onion states are separated by a vortex state. 
This is the dominant switching process for most re- 
ported magnetic rings, stabilized by any spurious 
spatial asymmetry. At the first transition point, 
cither one wall is pinned stronger than the other or 
both walls move towards each other due to asym- 
metry. In both cases the two domain walls ap- 
proach each other with increasing applied magnetic 
field. When these two domain walls meet, they an- 
nihilate each other, provided the ring is sufficiently 
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twisted states have recently been found to exist for 
small narrow rings^i. In most cases, the remanence 
state is an onion state. However, when the ring 
is sufficiently thick and wid o 29 i 30 , the remanance 
state is a vortex state. 

The third process is triple switching O-V-VC-O, 
shown in Fig|^d) i 30 i 31 ' 32 The switching mechanism 
is essentially the same as for O-V-O, except that 
the vortex state does not deform into an onion state 
abruptly, but first nucleates a vortex core in one 
arm of the ring. The core then moves slowly to 
the outer rim. During this part of the process the 
state is called a "vortex core" state (VC). To some 
extent, this is similar to the switching process of 
magnetic discs, where the core in the vortex state 
is shifted by the magnetic field and ultimately exits 
at the boundary (see Fig. Ufa)). 



II. SWITCHING PHASE DIAGRAMS 

We begin by summarizing our findings along with pre- 
viously reported results in the switching phase diagrams 
shown in Fig[TJ Until now, there are three dominant mag- 



These magnetization reversal processes are found in 
relatively large rings (Rq 3> L ex ) with spatial dimen- 
sions in the micrometer range. However, it is reasonable 
to expect that in the sub-micron and truly nano ranges, 
new phenomena exist and wait to be explored. There- 
fore, we have performed simulations for a wide range of 
parameters and identified the following three additional 
magnetization reversal processes. 
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FIG. 2: (color online) Magnetic switching processes in nanor- 
ing structures and their associated hysteretic behavior, (a) 
hysteresis of magnetic discs (Ri — > 0), (b) one-step O-O 
switching of rings, (c) double switching O-V-O of rings, (d) 
triple switching O-V-VC-O of rings. 



FIG. 3: (color online) Typical hysteresis curve for out-of-plane 
O v -V-0 switching (R — 70 nm, w — 7 nm, h — 10 nm). 
The insets are the snapshot during the O v -V transition. The 
areas with strong out-of-plane components are highlighted by 
open red circles. The big blue arrow points to a small step 
indicating the appearance of out-of-plane domains. 



B. Vertical Onion - Vortex - Onion Switching: 

o v -v-o 



On first sight, this process appears to fall in the 
same category as O-V-O, since its hysteresis, shown in 
Fig. [31 shares characteristics similar to double switching 
(Figj^c)). The switching mechanism behind it, however, 
is quite different. The onion state here does not consist 
of two transverse head-to-head domain walls, but rather 
of two vertical head-to-head domain walls (see the areas 
highlighted by open red circles in Fig. [3J , where the spins 
point out-of-plane rather than in-plane. Hence, we call 
the state involving vertical head-to-head domain walls a 
"vertical" onion state (O v ). The big blue arrow in Fig. 
[3j points to a small step in the magnetization curve, indi- 
cating the appearance of out-of-plane domains. The two 
domains may point parallel or antiparallel. However, the 
parallel configuration is preferred, since the switching of 
O-V is more controllable and the transition is sharper 
compared with the antiparallel case. As we can see from 
Fig. [3J when the two domains point in the same di- 
rection, the O-V transition is simply achieved by rotat- 
ing spins out-of-plane in the intermediate region between 
the two domain walls (highlighted by the red rectangle), 
which subsequently align with the out-of-plane domain 
wall spins. In the antiparallel case, the domain wall which 
moves faster will first rotate to form a normal in-plane 
head-to-head domain wall. Then the other domain wall 
rotates to in-plane alignment, such that the two domain 
walls can easily annihilate. This second process results 
in a slightly larger and broader transition field, which 
however can be easily avoided by applying a small out- 
of-plane magnetic field. 

The O v -V-0 process has not been reported before, as 
most of the initial attention in the field had been focused 



on thin film rings (h <C w). As seen in the phase dia- 
grams of Fig. [TJ this switching process indeed belongs 
to rings with w < h (see the yellow regions). Shape 
anisotropics force the spins in the domain walls to point 
out-of-plane. This process is crucial when one wants to 
use vortex states in high-density memory applications. In 
the thin film case (h < L ex ), vortex states are not easy 
to form when the ring radius is small (say R a < lOOnm). 
Instead, twisted states will stand in the way, preventing 
the formation of vortex states, which will be discussed 
in detail next. The only limit of this O v -V-0 process 
is the single domain limit for the ring. As long as onion 
states exist, this scheme applies. However, when the ring 
is extremely small, the onion state will be replaced by 
a single domain state. Roughly speaking, this transition 
happens at R ~ L ex . We use the remanence magne- 
tization as an order parameter to locate this transition, 
as shown in Fig. |4j The onion state is characterized by 
a relatively small magnetization (approximately 0.68 for 
the chosen w and h), since most of the spins are aligned 
along the ring boundary, whereas the single domain state 
is characterized by a magnetization close to saturation, 
since in this case nearly all the spins align in the same 
direction. Obviously, the transition region (shaded area 
in Fig. 0} is R a ps 20 - 40 nm w 1 - 2L ex . For rings 
with R Q < L ex , the only magnetization reversal process 
will be one-step switching via coherent single domain ro- 
tation. Therefore, the R of rings has to be at least L ex 
if the vortex state is desired. 
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FIG. 4: (color online) Single domain to onion state transition 
with incrreasing ring diameter. Here, the height is fixed at 
23 nm, and the width is fixed at 15 nm. The shaded area 
represents the transition region. 
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FIG. 5: (color online) Typical hysteresis curve for quadruple 
O-T-V-O switching (R = 150 nm,™ = 60 nm, h = 15 nm). 
The insets (a) and (b) are snapshots during the transition, (a) 
is the onion state at remanence, while (b) is a typical twisted 
state, characterized by a 360° domain wall. 



C. Quadruple Switching: O-T-V-O 

This switching process consists of four steps. In con- 
trast to regular O-V-O triple switching, an additional 
step, characterized by the appearence of a twisted state 
(see Fig. [5] inset (b)), is involved. As mentioned above, 
twisted states are only found in rings with relatively small 
lateral dimensions 2 ^ (R < 500nm). In order to under- 
stand the conditions for the twisted state to occur, first 
recall the process of onion-to-vortex (O-V) transition 
during the regular triple switching: two domain walls 
move towards each other and annihilate by introducing 
a vortex core passing through the rim, i.e. entering from 
the outer edge and exiting through the inner edge. As 



the diameter of the vortex core is typically equal to L ex , 
the width of the ring needs to be at least of the order of 
several exchange lengths in order to accommodate such 
a process. However, when the ring width is sufficiently 
small, the annihilation process is hindered. In that case, 
the two domain walls are stuck together, forming a 360° 
domain wall (see Fig. [5] inset (b)). The configuration in- 
volving this 360° domain wall is called "twisted state" i^i 
It turns out that this state can be stabilized within a 
fairly large field range AH, e.g. for the specific ring 
shown in Fig. [5] it is stable between 100 Oe and 300 
Oe. This result is consistent with previous literatur o 21 ' 23 . 
Moveover, the results in Ref. Hi] show precursor behav- 
ior of this type of switching, since the vortex state will 
not give positive magnetization when the magnetic field 
is negative. 

From our simulations we find that the field range AH 
of the twisted state is sensitive to the geometric param- 
eters (i? , w, h) in the following ways. 

1. h < L ex : We find that the twisted state is not 
stable when the ring height is larger than L ex . It 
has been reported that magnetic domain structures 
vary significantly along the vertical direction when 
the height is l&rgeM- Our simulations show that 
this is indeed the case when h > L ex , as the ex- 
change length L ex is a measure of typical domain 
size. In this case, the vertical stacking of slightly 
off-set domains eases the introduction of vortices in 
the annihilation process, destabilizing the 360° do- 
main wall. Furthermore, AH is most sensitive to 
ring height. The flatter the system, the larger AH, 
i.e. AH can be doubled with decreased height. On 
the other hand, vortices cannot be formed if h is too 
small (< Q.3L ex ), which leads to the last switching 
process we will discuss in the next subsection. 

2. L ex < w < 5L ex : Typically, the field range AH 
varies with ring width (about 20%), with a peak 
around w = 3L ex . The reason for this is that 
there is a qualitative difference between thin and 
wide rings when the 360° domain wall is destroyed. 
For rings with w < L ex , the 360° domain wall is 
not stable, with the exception of ultra-small rings 
(R < 2L ex ). As it is the case for out-of-plane O- 
V-O switching, the spins in the domain area first 
point out-of-plane, and then turn in-plane along the 
field with increasing applied field. Slightly wider 
rings can support the 360° domain wall in a small, 
but finite, field range and then evolve with the 
same process. For wider rings with w > 5L ex , a 
vortex can be introduced through the outer edge 
of the ring and traverse through the rim, as al- 
ready described for the O-V-O triple switching 
case. Slightly thinner rings still feature vortices. 
However, in this case the vortices cannot remain 
in the area of the ring where they entered when 
the reversed transverse field is increased. Instead, 
they move along the rim of the ring and finally 
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exit. Hence, in this case the associated magnetiza- 
tion step in the hysteresis curve has a fairly large 
slope, indicating that the twisted region is moving 
with increasing applied field. The competition of 
these two mechanism gives rise to a peak of the 
field range where the twisted state is stable, which 
occurs at w ~ 3L ex . 

3. L ex < R a < 15L ex : In order to observe the twisted 
state, the entire ring has to be sufficiently small, 
but still larger than the single domain limit. The 
reason why larger rings destabilize the twisted state 
more easily is that in these structures vortices can 
enter from the outer rim. In smaller rings, vor- 
tices are not as easily accomodated in the domain 
wall region, because of differences between the ar- 
clengths of the outer and inner domain wall bound- 
aries, due to the increased curvature. Another rea- 
son is that for large rings, there is a higher chance 
that defects serve as a catalyst of local vortex for- 
mation. AH is almost independent of R for fixed 
aspect ratios (k = Ri/R a ). On the other hand, the 
width of the magnetization step in M(H) where the 
vortex state is stable decreases quickly when R a is 
decreased, because the curvature eases the process 
of introducing a vortex during the V-O process. 
When k is large and the ring is ultra-small, the ex- 
treme case where no vortex is involved emerges as 
discussed in detail below. 

The above conditions are fulfilled in the gray regions 
of the phase diagrams in Fig. [TJ 




FIG. 6: (color online) Twisted triple switching (O-T-O): (a) 
Hysteresis curves for ring 1 (_R = 70 nm, w = 18 nm, h — 5 
nm, black continuous line) and ring 2 (R = 42 nm, w — 16 
nm,h — 14 nm, red dashed line); (b) spin configuration for 
ring 1 at H — 200 Oe; (c) spin configuration for ring 1 at 
H = 600 Oe; (d) spin configuration for ring 1 at H — 650 Oe. 



D. Twisted triple switching (O-T-O) 

This process is an extreme limit of the O-T-V-O 
switching process. Rather than forming a vortex state to 
lower the energy, the system chooses to move the entire 
domain wall towards either end of the ring along the mag- 
netic field (see Fig. (He)), when the applied magnetic field 
is increased. Ultimately, the 360° domain wall breaks 
up into two 180° domain walls, and one of them travels 
back to the other side (see Fig. [6jd)), forming the re- 
versed onion state. For relatively large rings, one can ob- 
serve a small step in the magnetization curve, indicating 
the occurence of the twisted state (see Fig. HJb)), while 
this step is not noticeable for very small rings. After 
the twisted state is formed, the magnetization decreases 
approximately linearly with increasing applied magnetic 
field, indicating the movement of the domain wall (see 
Fig. [UJc)). Finally, a large jump is observed, indicating 
the breakup (see Fig. [S^d)) of the 360° domain wall and 
the formation of the reversed onion state. This process 
can only be observed under two extreme conditions: 

1. Very small narrow rings (typically R < 3L ex and 
K > 0.5), when 0.3L ex < h < L ex (see the small 
red region in Fig. HJa)). The red dashed line in 
FigEJa) is the typical hysteresis in this category. 
Here it is the curvature of the ring that prohibits 
the annihilation process. Under such conditions, 
the domain region typically covers a quarter of the 
ring. 

2. Ultra-thin small narrow rings (typically h < 0.3L ex 
and k > 0.5). The black solid line in FigH^a) is 
the typical hysteresis. Fig[6jb)(c)(d) are snapshots 
of the configuration evolution with increasing ap- 
plied magnetic field. In this case, it is the strong 
easy-plane shape anisotropy that hinders the an- 
nihilation process. For these ultra-thin systems, 
this switching process can be observed for rings as 
large as R ~ lOOnm (see the fairly large red region 
in FigQTc)). Theoretically they may exist on even 
larger rings. In reality, however, local defects will 
favor the formation of a local vortices, destroying 
the 360° domain wall. 

Note that the twisted state is stable at remanence once 
it is formed. Designs utilizing this state for magnetic 
memory have been proposed^. Our results show that 
such designs can use very small rings, for which there 
appears to be great potential to build high density stor- 
age systems. 

III. CONCLUSIONS 

In conclusion, we have computed the phase diagrams 
of magnetic nanorings, and identified three new types of 
magnetization reversal processes, such that there are now 
six known types of magnetization reversal processes in all 
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reported for magnetic nanorings. The O v -V-0 switch- 
ing features a new out-of-planc onion state and easy 
switching process. Rings in this category have h > w, 
and R Q can be as small as L ex . The quadruple switch- 
ing (O-T-V-O) and twisted triple switching (O-T-O) 
feature twisted states. Rings in these two categories are 
fairly small (typically with outer radius less than 200nm) . 
Their width has to be less than about AL ex to stabilize 
the 360° domain wall. The field range where the twisted 
state is stable is most sensitive to the ring height. Small 
heights give rise to large field ranges, and there will be 
no twisted state when h > L ex . We have combined our 
findings with results from previous literature into com- 
plete switching phase diagrams, shown in FigQ] These 
should be helpful in designing magnetic nanorings with 



specific properties, and show the potential of building 
high density storage systems based on ring structures. 
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